Proton and pion transverse spectra at RHIC from radial flow and finite size eff'ects 
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We show that the proton and pion transverse momentum distributions measured at RHIC for all 
collision centralities for pions and most of the collision centralities for protons, can be simultaneously 
described in terms of a thermal model with common values for the radial flow and temperature, 
when accounting for the finite size of the interaction region at the time of decoupling. We show that 
this description is obtained in terms of a simple scaling law of the size of the interaction region with 
the number of participants in the collision. The behavior of the proton to pion ratio at mid-rapidity 
can also be understood as a consequence of the strength of the radial flow and system size reached 
at RHIC energies. 



PACS numbers; 25.75.-q 



I. INTRODUCTION 



The unexpected behavior of the proton to pion ratio as 
a function of has been taken as an indication of the on- 
set of the thermal recombination of quarks as an impor- 
tant mechanism for hadron production at RHIC energies 
in Au + Au collisions In its simplest form, thermal 
recombination invokes a densely populated parton phase 
space to allow the statistical formation of hadrons from 
constituent quarks assigning degeneracy factors appro- 
priate for either mesons or baryons 

An important short come of the recombination sce- 
nario is that it ignores inelastic and elastic scattering ex- 
perienced by hadrons before kinetic freeze-out and thus 
that neither particle abundance nor their spectra are 
fixed right after recombination. A more appropriate de- 
scription of statistical systems, that includes the fact that 
the detailed history is washed out by means of interac- 
tions after a long enough time, can be given in terms of 
global features which survive all the way to the end of 
the system's evolution. One of these global features is 
fiow, in particular radial flow. 

It has been observed that the magnitude of the radial 
fiow velocity exhibits a 50% increase from AGS and SPS 
to RHIC energies Recall that in p+p collisions, where 
no effects of radial flow exist, it is known that the proton 
to pion ratio as a function of pt remains basically un- 
changed, never exceeding one, for collision energies rang- 
ing from 19.4 GeV at the Tevatron, 44.6 and 52.8 GeV at 
ISR up to 200 GeV at RHIC 4]. In contrast, the proton 
to pion ratio in Au + Au collisions at RHIC reaches and 
even exceeds one for pt ^ 2 GeV. Therefore, if recombi- 
nation of thermal partons has anything to do with this 
behavior it is clear that a thermal description of the in- 
dividual particle spectra has to be possible at least up to 
such Pt values. Nevertheless, in Ref. , a fit to a thermal 
model that attempts a description of particle spectra up 
to Pt ^ 2 GeV in terms of an intrinsic freeze-out temper- 
ature To, together with radial flow, yields values of order 



Tq ~ 180 MeV, which is closer to the hadronization tem- 
perature than to the kinetic freeze-out temperature. 

Yet another intriguing behavior that concerns freeze- 
out temperatures and expansion velocities in thermal 
models is their relation as a function of the centrality 
of the collisions, which, within the usual thermal model 
calculations, can be stated as an increase in fiow together 
with a decrease in temperature as the centrality of the 
collision increases. This behavior is usually attributed 
to the larger time spent by the system in the hadronic 
phase for the most central collisions allowing for the de- 
velopment of fiow and consequently decreasing the values 
for the kinetic freeze-out temperature Q • However, as 
results from elliptic fiow analyses seem to indicate [1], 
fiow is generated early, in the partonic phase of the colli- 
sion. Moreover, kinetic freeze-out temperatures can also 
be thought of as a global feature of strongly interacting 
systems that refiect the average kinetic energy needed for 
the system to decouple. In fact, a systematic study of 
HBT data and particle yields for pions at mid-rapidity 
from AGS to RHIC energies Q shows that this aver- 
age energy is independent of centrality and beam energy. 
Therefore one can ask if an alternative description, with 
common values of temperature and flow velocity, reflect- 
ing the above property of strongly interacting systems, 
can be achieved for all centralities. As we show, the key 
ingredient that allows a description of particle spectra 
in a thermal model, including radial fiow, and that ad- 
dresses the above mentioned phenomena from a unifying 
point of view, is the realization that particle production 
and successive freeze-out in a relativistic heavy-ion envi- 
ronment takes place during small time scales, of order 10 
fm, and consequently within small volumes. 

Although not commonly considered, small size effects 
are important in the description of a variety of phenom- 
ena associated with statistical systems such as the late- 
stage growth of nucleated bubbles during a first order 
phase transition jTd| and the statistical hadronization 
model Finite size effects are also known to infiu- 
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ence the interpretation of the correlation lengths in Han- 
bury Brown- Twiss anal ysis in the context of rclativistic 
heavy- ion collisions [l^. Il3l| . 

Recall that useful microscopical information in this 
kind of collisions can be obtained from comparing the 
average interparticle separation during the collision evo- 
lution to the range of strong interactions. For the case 
of pions (the most copiously produced particles in the 
collision) right after the collision, the system is better 
described as a liquid rather than as a gas One im- 
portant consequence is the appearance of a surface ten- 
sion which acts as a reflecting boundary for the parti- 
cles that move toward it. The reflection details depend 
on the wave length of the incident particle but the im- 
portant property introduced by the reflecting surface is 
that it allows very little wave function leakage and, to a 
good approximation, the wave functions vanish outside 
this boundary. When the average separation of the par- 
ticles in the system becomes larger than the range of the 
strong interaction, they become a free gas but because 
of the short interaction range, the transition between the 
liquid and the gas stages is very rapid and the momen- 
tum distribution is determined by the distribution just 
before freeze-out. 

To be concrete, we need to compare the pion separation 
d to the average range of the pion strong interaction (ds ~ 
1.4 fm). For typically accepted values for the density and 
formation times @, it is possible to show that d ^ 0.6 
fm < ds and the condition to regard the pion system as 
a liquid is met. 

Qualitatively, the behavior of thermal particle spectra 
including finite size effects deviates from a simple expo- 
nential fall-off at high momentum since from the Heisen- 
berg uncertainty principle, the more localized the states 
are in coordinate space, the wider their spread will be 
in momentum space. In terms of the discrete set of en- 
ergy states describing the particle system, this behavior 
can be understood as arising from a higher density of 
states at large energy as compared to a calculation with- 
out finite size effects. These ideas have been applied to 
the description of charged and neutral pion spectra mea- 
sured at RHIC with a good agreement for the transverse 
momentum interval < pt ^ 3 (isj . 

In this paper we compute the transverse momentum 
distribution for pions and extend the above ideas to 
also include protons in the description, assuming ther- 
mal equilibrium together with radial fiow and account- 
ing for finite size effects at decoupling. By comparing to 
data on pion and proton spectra on Au+Au collisions at 
^SNN = 200 GeV y , we show that for temperatures and 
collective transverse flow within values corresponding to 
kinetic freeze-out conditions, the transverse momentum 
distributions can be described with common values of 
temperature and expansion velocities for all collision cen- 
tralities for pions and for most of the collision ccntralities 
for protons. 

The work is organized as follows: In Sec.|TT]we present 
the basics of the model to compute pion and proton dis- 



tributions. In Sec. mil we compute the transverse mo- 
mentum distributions for protons and pions comparing 
to data on Au+Au collisions at y^s^N = 200 GeV. We 
show that a good agreement with these data for different 
collision centralities can be achieved by assuming a sim- 
ple scaling of the radii with the cube root of the number 
of participants in the collision. In Sec. IIVI we compute 
the pion correlation function and also extract the size of 
the system as a function of the cube root of the number 
of participants in the collision. We flnally conclude in 
Sec. El 



II. THE MODEL 

We consider a scenario where flnite size effects are in- 
cluded by restricting the system of particles to be con- 
flned within a volume of the size of the flreball at freeze- 
out. Since we aim to describe spectra at central rapidi- 
ties, it suffices to take the confining volume as a sphere 
of radius R (fireball) as viewed from the center of mass of 
the colliding nuclei at the time of decoupling [16,]. This 
time needs not be the same over the entire reaction vol- 
ume. Nevertheless, in the spirit of the fireball model we 
consider that decoupling takes place over a constant time 
surface in space-time. This assumption should be essen- 
tially correct if the freeze-out interval is short compared 
to the system's life time. Though some of the parti- 
cles emitted in the central rapidity region could originate 
from a finite range of longitudinal positions due to ther- 
mal smearing, we will consider that most of the central 
rapidity particles come from the central spatial region 
and thus neglect possible effects on these particles from 
a different longitudinal and transverse expansion veloci- 
ties. 

In the case of bosons, the wave functions that incor- 
porate the effects of a finite size system have been found 
in Ref. 16], where we refer the reader for further details 
of the model. These wave functions are given as the sta- 
tionary solutions of the wave equation for bosons, namely 
the Klein- Gordon equation 
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(r,i) =0 



(1) 



subject to the boundary condition (j>{\r\ — R,t) — , and 
finite at the origin. The normalized stationary states are 



(r, t) 
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In the case of fermions, the wave functions are found as 
the stationary solutions of the Dirac equation 



(i7^a^-m)V(r,t) = 



(3) 



subject to the the boundary condition V'(|r| = R,t) = , 
and also finite at the origin. It is easy to show that the 
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normalized stationary states are 



2mR Ji+3/2ikniR) 
Eni + m + i<T ■ 

— Enl + m - Z(T • V 

Yi„i'{r)Ji+i/2{knir) 



(4) 



In Eqs. ([2]) and ([4]), Ji, is a Bessel function of the first 
kind and Yi„i' is a splierical harmonic, cr are the Pauli 
matrices and the parameters kni are related to the energy 



eigenvalues E^i by 



and are given as the 



solutions to Ji+i/2{kniR) — 0. The contribution to the 
thermal invariant distribution from a state with quantum 
numbers {n, I, m'} is given by 



E 



(2^ 



{knl ■ u)f{knl ■ v)yVnlm'{P,r), (5) 



where Wnim' (p, r) is the Wigner transform and /(fc„; • v) 
the thermal occupation factor of the state, respectively. 
The four-vectors u'' and represent the collective flow 
four-velocity and a four-vector of magnitude one, normal 
to the freeze-out hypersurface S, respectively. 

In order to consider a situation where freeze-out hap- 
pens at a fixed time and within a spherical volume of 
radius R, the unit four-vector u'^ can be chosen as u'^ = 
(1,0). To keep matters simple, we also consider a ther- 
mal occupation factor of the Maxwell-Boltzmann kind 
fikni • v) — e^*^"' "/^ where T is the system's tempera- 
ture. The four-vector is parametrized as v'^ = 7(1, v), 
and we choose a radial profile for the vector v such as 
V — (3r/R, where the parameter (3 represents the surface 
expansion velocity. Correspondingly, the gamma factor 
is given by 



1 



(6) 



Nonetheless, in order to continue to keep matters as 
simple as possible, and be able to analytically perform 
the integrations in Eq. ([5]), we will instead consider that 
the gamma factor is a constant evaluated at the average 
transverse expansion velocity, namely, 



7^7 = 



v/1- (3/3/4)= 



(7) 



where the average is computed by assuming that the mat- 
ter distribution is uniform within the fireball. 

We take the four- vector fc^^ — k„;), and choose 
k„i |j p. This choice is motivated from the contin- 
uum, boundless limit, where the relativistically invari- 
ant exponent in the thermal occupation factor becomes 
7(£;-p-v). 
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FIG. 1: (Color online) Invariant tt"*" and p distributions as a 
function of pt for T = 117 MeV, i? = 8 fm and = 0.6 
and /3p = 0.53. Also shown in the insert is the ratio pj-K^ of 
these distributions. Data are from Ref. [J] for collisions with 
- 10% of centrality. 



Summing over all the states, the invariant thermal dis- 
tribution for bosons and fermions are given by 
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respectively. The factor {21 + 1) in Eqs. ^ and ([9]) comes 
from the degeneracy of a state with a given angular mo- 
mentum eigenvalue I. A/jj and A/f are normalization con- 
stants. 
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FIG. 2: (Color online) Invariant tt^ (left panel) and p (right 
panel) distributions as a function of pt for different collision 
centralities. The uppermost curves correspond to the most 
central collisions, — 10% with the centrality decreasing in 
intervals of 10% Jj] multiplied by successive factors of 10~^. 
The size of the equivalent spherical region is calculated ac- 
cording to the relation R = Ro + C{Np^rt/2y^^ , wfth Ro = 1 
fm and C = 1.28. The curves represent the theoretical calcu- 
lation with Tp,n = 117 MeV, Pn = 0.6 and /3p = 0.53, which 
describe the most central collisions data. 



The contrast between a calculation with and without 
finite size effects, can be appreciated by looking at Fig. 5 
in Ref. [isj were it is shown a comparison between the 
invariant pion distribution as a function of pt computed 
for T = 120 MeV, /? = 0.6 with finite size effects (i? = 8 
fm) and without them. The curves are also compared 
to data on positive pions from PHENIX [l|. We notice 
from the figure that the curve with finite size effects does 
a very good job describing the data for all values of pt in 
this range. In contrast, a calculation where no effects of 
a finite size are included, and thus the wave function of a 
given state is simply a plane wave, does not describe the 
data over the considered range when use is made of the 
same values for T and /3 as for the case of the calculation 
with a confining volume. 



III. TRANSVERSE SPECTRA 

We now compare the model to data on mid-rapidity 
positive pions together with protons from central Au + 
Au collisions at JJnn = 200 GeV measured in RHIC Q. 
We perform a x fit to each spectra. The fit parameters 
are the pion and proton fireball radii, R-j^, Rp, temper- 
atures Tjr, Tp, surface radial flow velocities /Syr, /3p and 
normalizations A/'tt, Np. Based on the success of the de- 



centrality 




R (fm) 


- 10% 


325.2 


8.0 


10 - 20% 


234.6 


7.3 


20 - 30% 


166.6 


6.6 


30 - 40% 


114.2 


5.9 


40 - 50% 


74.4 


5.3 


50 - 60% 


45.5 


4.6 


60 - 70% 


25.7 


4.0 


70 - 80% 


13.4 


3.4 


80 - 92% 


6.3 


2.9 



TABLE I: Parameters A'part and R for the description of the 
pion and proton spectra corresponding to different central- 
ities. The radii of the equivalent spherical region has been 
calculated according to the relation R = Ro + C{Np.^^t 12)^^^ , 
wfth i?o = 1 fm and C = 1.28. 



scription of the central rapidity pion data obtained in 
Ref. [Tsil up to Pt ~ 3 GeV, we first fix the parameters 
describing the pion data with the minimization proce- 
dure. The parameters thus obtained are R-j^ — 8 fm, 
Tyr — 117 MeV and jS-u — 0.6 which are basically the 
same as the ones obtained in Ref. [isj were only the nor- 
malization was left as a free parameter and the rest were 
set to reasonable values that describe freeze-out condi- 
tions at RHIC. 

Next, in order to find the parameters that describe 
the proton spectrum, we fix the values of any two of 
the parameters Rp, Tp and /3p to be the same as the 
corresponding parameters describing the pions, leaving 
the third parameter, along with the normalization con- 
stant A/p free. The optimum set of parameters obtained 
with this procedure correspond to Rp = 8 fm, Tp = 117 
MeV and /3p = 0.53. Figure [1] shows the distributions 
for pions and protons for central collisions (0 — 10%) [l| 
compared to the theoretical calculation with the best pa- 
rameters obtained. We notice that the proton data are 
well described by the model up to pt ~ 3 GeV for a tem- 
perature and system's size equal to the corresponding 
parameters for the pions but that the magnitude of (3p 
is about 10% smaller than We recall that in order 
to find an analytical expression for the momentum dis- 
tributions, we resorted to approximate the gamma factor 
in Eq. ([6]) by the average gamma factor in Eq. ([7]) . Since 
the effect of the same radial fiow is stronger for parti- 
cles with larger mass, it is therefore natural to expect 
that with this approximation we introduce a discrepancy 
in the description of the ffow for particles with different 
masses. Nevertheless, we feel that an error of order 10% 
is acceptable considering the advantage of working with 
analytical expressions. 

We also notice that the description of the proton data 
for Pt > 3 GeV is not as good. This can be understood 
recalling that for large pt the leading particle production 
mechanism is the fragmentation of fast moving partons, 
some of which fragment outside the fireball region and 
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thus are not influenced by the confining boundary that 
the rest of the particles experience within the fireball and 
thus, that our description is not valid for these large pt 
particles. 

Figure [2] shows the distributions for pions and protons 
for difi'erent collision centralities. For the description of 
these data, we have fixed the values of Tp^ ti and (3p^ -j^ to 
the ones obtained from the most central collisions anal- 
ysis leaving the normalizations to be determined by the 
minimization procedure. The size of the overlap region 
for peripheral collisions has been determined from the 
number of participants N^^^t in the reaction [H by a sim- 
ple scaling law for the size of the equivalent spherical 
region according to the relation i? = i?o + C'(-^part/2)^/"^, 
with i?o = 1 fm and C — 1.28 that gives i? = 8 fm for 
the most central region (0 — 10%) data. This relation is 
motivated by the similar one that gives the radius of a 
nucleus in terms of the mass number. The value of i?o 
tries to account for the finite size of the interaction region 
as the number of participants takes its smallest value for 
the most peripheral collision, namely, 7Vp„t = 2. 

The values for R and A^part are listed on Table HI We 
notice that the pion data are well described for all cen- 
tralities except at the lower end of the spectra where a 
thermal calculation is expected to fail due to resonance 
contamination. The proton data is well described up to 
Pt ~ 3 GeV only up to centralities of order 40 — 50% from 
where the quality of the description decreases as the cen- 
trality of the collisions decreases. We interpret the poor 
description of the proton data for pt > 3 GeV for all cen- 
tralities as an indication that the leading particle produc- 
tion is not a kind of thermal parton recombination but 
instead the fragmentation of fast moving partons. On 
the other hand, the failure to describe proton data for 
centralities smaller that 40 — 50% could be attributed to 
a different scaling of the effective size of the interaction 
region with iVp^,! as compared to the one obeyed by the 
pions or to the fact that for large impact parameters, the 
proton size becomes comparable to the size of the inter- 
action region. An analysis to explore these possibilities 
will be presented elsewhere. 



IV. PARTICLE CORRELATIONS 

In order to explore the space-time dimensions of the 
system created in high-energy heavy-ion collisions, one 
typically looks at two-particle correlation functions. In 
the present case, it is thus instructive to look at this func- 
tion to see whether the size scale that can be extracted 
from a two-particle correlation is comparable to the in- 
trinsic scale dimension used in the model formulation. 
We carry out the analysis for the two-pion correlation 
function. For the purposes of this section, we closely fol- 
low Ref. [l2] to where we refer the reader for details. 

Let tpnim'ip) represent the Fourier transform of the 
wave function for the state with quantum numbers 



n, l.m' , namely 



tpnlm' (P) 



-e-*P-^Vnim'(r). 



(10) 



With the normalization adopted in Eq. ([2]), the one-pion 
momentum distribution can be represented as 



A(P) 



d^p 

2£;„, 



te 



Cim'(P)V'nim'(p) -(11) 



n,l,m' 



Similarly, and under the assumption of a complete factor- 
ization of the two-particle density matrix, the two pion 
momentum distribution can be written as 



-P2(Pl,P2) 



d^N 



d^pid^P2 

Pl{Y>l)Pl{P2) 



m' (P2 ) 



nlm' 



(12) 



from where the two-pion correlation function C2 can be 
written, in terms of Pi and P2 , as 



C2(Pl,P2) 



f2(Pl,P2) 
A(Pl)A(P2)' 



(13) 



Notice that as a consequence of the factorization assump- 
tion, the correlation function is such that C2(p,p) — 2. 

For the spherically symmetric problem here described, 
the correlation function depends on the magnitude, as 
well as on the angle between the momenta of the two 
particles pi and P2. We make the change of variables to 
relative q = Pi — P2 and average K = (pi -I- p2)/2 mo- 
menta and also to the angle between these last two vec- 
tors, 9. The correlation function thus becomes a function 
ofiir=|K|,(7=|q| and 9. In order to consider the con- 
tribution from pions with different angles between their 
momenta, we average over 9. Figure 3 shows C2{q) av- 
eraged over 9 and for a fixed value K — 260 MeV as a 
function of g for i? = 8 fm, T = 117 MeV and 13 = 0.55. 
The solid curve shows the corresponding Gaussian fit. 

In order to extract the system's size R^s from this func- 
tion, we fit this curve to a Gaussian distribution of the 
form 



5(9) 



1 + P'('Z), 
exp(-(72i?2^/2). 



(14) 



Figure 4 shows the behavior of i?„ff as a function of 
(A^part)^^'^ compared to measured values for 7?^;^^ [l3] from 
Au + Au collisions at ^Jsmn = 200 GeV. The lower solid 
curve corresponds to our model effective radii whereas 
the upper solid curve is the model curve displaced by a 
constant R = 0.8 fm. The dashed curve corresponds to 
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FIG. 3: (Color online) C2(g) averaged over the angle between 
K and q as a function of q for K — 260 MeV, 7? = 8 fm, 
T — 117 MeV and /3 = 0.55. The histogram corresponds 
to the model calculation whereas the solid curve shows the 
corresponding Gaussian fit. 

the best linear fit to the data. We notice that the slope 
of our model curve is in good agreement with the data. 
The fact that the intercept is different from zero may in- 
dicate the existence of correlation scales in the data that 
are not considered in our simple approach. 
Finally, recall that 

p{r) - exp(-rV2i?L) (15) 

is the spherically symmetric three-dimensional distribu- 
tion in space that gives rise to p{q) upon Fourier trans- 
formation and that the Rr m s radius from p{r) is given 

by 

Rr.rn.s. = ViRoil- (16) 

On the other hand, for a rigid sphere, such as the distri- 
bution giving rise to our model distribution, the Rr.m.s. 
radius is given by 

Rr.r,^.s. = V^JbR. (17) 

By equating these two r.m.s. radii, we see that in order 
to compare the effective radius with the model one, the 
relation between them is given by 

R = VSEeff (18) 



V. CONCLUSIONS 

In conclusion, we have shown that by considering the fi- 
nite size of the interaction region and a simple scaling law 
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FIG. 4: (Color online) R^a as a function of (A^p„t)^''^ for 
K — 260 MeV compared to measured values for iiaido from 
An + Au collisions at y^sjvjv ~ 200 GeV. The lower solid 
curve corresponds to the model effective radii whereas the 
upper solid curve is the model curve displaced by a constant 
Ro ~ 0.8 fm. The dashed curve is the best linear fit to the 
data 



of this with the number of participants in the collision, it 
is possible to achieve a good description of mid-rapidity 
pion and proton data in Au -I- Au collisions at RHIC with 
common values of temperature and transverse expansion 
velocity. For central collisions, the proton to pion ratio 
is also well described and its behavior can be attributed 
to the strength of the radial flow achieved in RHIC. By 
performing a two-particle correlation analysis and com- 
paring to data for R^^^ia as a function of (iVp^t)"'^^'^ from 
Au -I- Au collisions at ^/snn — 200 GeV, we see that the 
scaling law found from the single particle spectra analysis 
is in good agreement provided that we displace our model 
curve by a constant R = 0.8 fm and we speculate that 
this signals the existence of an extra correlation length 
in data that is not accounted for in our simple model. 

We should stress that the spherical symmetry assumed 
throughout can be thought of as a theoretical tool rather 
than as a realistic approximation to the actual collision 
geometry for the highest RHIC energies. Our intention is 
to provide a working model with a high degree of symme- 
try that can be better controlled in an actual calculation. 
The same is true for the treatment of the shape of the 
region for non-central collisions which one knows, from 
elliptic flow analyses, that partially retains the original 
almond shape of the overlap region in the collision. The 
sizes referred to in this way reflect characteristic sizes 
rather than actual spherical radii. Although modifying 
the geometry used for the calculation will certainly give 
rise to a different set of quantum states, the bulk of the 
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effect will remain since the physics that it captures is the 
Heisenberg uncertainty principle whereby restricting the 
size of the region to become finite, the momentum states 
become broader. As for the use of stationary states, we 
point out that although the freeze-out volume is reached 
with a large expansion velocity, the transition from a 
strongly interacting system to a free gas is rapid and 
what matters is the distribution right before this transi- 
tion and therefore the length scale associated to it. 

It is important to emphasize that a description of the 
transverse distributions for different impact parameters 
can be done by considering a varying freeze-out temper- 
ature and radial velocity but the lesson to be learned 
from the present analysis is that this variation can be 
tempered and/or even avoided by considering the finite 
size of the interaction region. 

We point out that some hydrodynamical models with- 
out finite size effects have been able to give similarly good 
descriptions of data up to pt of order 2-2.5 GeV jjj] at 



the expense of introducing a large amount of parame- 
ters. What we have shown here is that it is also possible 
to achieve the same quality of description including a ba- 
sic property of quantum systems often neglected, that is 
the fact that in high energy reactions, particles are pro- 
duced in small space-time regions. While doing this and 
in this first step approach, we made use of approxima- 
tions to render the calculations tractable, one of such is 
the treatment of the gamma Lorentz factor in terms of an 
average one. The relaxation of these approximations is a 
natural step forward and we will report on the progress 
of this work elsewhere. 
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